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In the frame of this paper a wind climatology has been built for Greece, based on high resolution model 
simulations performed for a typical year of wind conditions over the area. The methodology followed 
includes: (a) the development of a typical wind year, (b) the definition of a modeling strategy and the 
performance of the appropriate simulations, (c) the verification of the 10-m wind speed simulated by the 
model against observations, (d) the statistical analysis of the wind simulations for the construction of 
the wind atlas. The analysis is performed at 50 m and includes gridded wind speeds, the parameters of 
the respective Weibull distributions, and the potential power production. Finally a discussion on the 
effect of the topography height and slope on the wind distribution and consequently on the wind 
potential density is provided. 
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1. Introduction 


Wind energy offers significant potential for near- and long-term 
carbon emissions reduction. The wind power capacity installed 
by the end of 2009 was capable of meeting roughly 1.8% of the 
worldwide electricity demand, and that contribution could grow to 
in excess of 20% by 2050 if reducing greenhouse gas emissions was 
to be considered as a high priority goal and considerable efforts 
were made to address impediments to increased wind energy 
deployment [1]. During the last decade the wind energy market 


* Corresponding author. Tel.: +30 210 8109126; fax: +30 210 8103236. 
E-mail address: kotroni@noa.gr (V. Kotroni). 


1364-0321/$-see front matter © 2013 Elsevier Ltd. All rights reserved. 
http://dx.doi.org/10.1016/j.rser.2013.10.016 


has expanded rapidly. From a cumulative capacity of 14 GW by the 
end of 1999, global installed wind power capacity increased 12-fold 
in 10 years to reach almost 160 GW by the end of 2009. Regionally, 
Europe continues to lead the market with 76 GW of cumulative 
installed wind power capacity by the end of 2009, representing 48% 
of the global total [2]. 

The Hellenic Government has set a target of 20% share of 
renewable energy in the gross final energy consumption by 2020. 
The contribution of renewable energy sources (RES) to the national 
energy balance in 2008 was approximately 7.8% of the gross final 
energy consumption and around 16.3% of the primary energy 
production. From the 1.64 Mtoe of primary energy produced from 
RES in 2008, 193 ktoe (~ 11.7%), have been produced by wind 
power plants [3]. The installed capacity of wind farms in Greece in 
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2011 was 1635.87 MW out of the 2511.60 MW of total installed 
capacity of renewable energy resources [4] while the installed 
capacity of wind power plants by 2020 is expected to exceed 
7000 MW. 

Based on the aforementioned information and taking into 
consideration the increase in the wind energy market, the identi- 
fication of potential sites for wind farms is of high interest and 
therefore the need of construction of a wind atlas over Greece is 
considered to be of high priority. The complexity of the terrain of 
Greece arises from its orography but also from the land-water 
distribution and the large number of islands. Therefore the 
production of an accurate wind atlas is also a quite complex task. 
Indeed, a wind atlas for Greece has been produced in 2006 [5] and 
updated in 2008, based on the use of measured winds over the 
Greek territory and the application of a mass-consistent model. 
Results of this wind atlas are given at the official web page of the 
Greek Regulatory Authority for Energy (http://www.rae.gr/geo/). 

Many wind atlases that have been recently produced for various 
countries have been based on the use of numerical weather predic- 
tion models (NWP) which are accompanied in some cases with the 
use of diagnostic downscaling models. A review of recent efforts in 
producing model based wind atlases, including the existing wind 
atlas for Greece, is given in Section 2 of this paper. 

Thus, based on the international practice to use NWP models 
for the construction of wind atlases, the National Observatory of 
Athens took the initiative to build a wind climatology for Greece 
based on high resolution model simulations performed for a 
typical year of wind conditions over Greece. This effort is con- 
sidered as complementary to the existing wind atlas that is based 
on the use of a mass-consistent model. The adopted methodology 
consists in defining a typical wind year over Greece that then is 
simulated using a mesoscale model for a domain covering Greece. 
In order to better represent the complex topography of Greece and 
the land-sea distribution a fine horizontal grid increment of 
2 km x2 km has been selected. The simulated wind fields have 
been verified against surface station observations in order to 
assess the model ability to reproduce the local wind flow. Finally 
the model simulated winds at various heights have been statisti- 
cally analyzed for the construction of the wind atlas. 

The rest of the paper is structured as follows. The next section 
presents a review of model based wind atlases presented in the 
literature, while in Section 3 the methodology adopted for the 
definition of the typical wind year over Greece is presented. 
Section 4 is devoted to the presentation of the modeling strategy 
and the verification of the model simulations, while Section 5 
discusses the resulting wind fields that are used for the construc- 
tion of the wind atlas. The last section of the paper is devoted to 
the concluding remarks of this work. 


2. Brief review of model based wind atlases. 


Many wind atlases that have been recently presented in the 
literature have been based on the use of either NWP models, or 
microscale models such as Computational Fluid Dynamics (CFD) 
models or mass-consistent type models, or a combination of 
mesoscale/microscale models. Indeed, the Swedish wind atlas 
has been produced based on the simulation of 192 samples 
selected to represent the meteorological conditions governing 
the wind climate of Sweden, using the CFD hydrostatic model of 
the Meteorological Institute Uppsala University (MIUU). The 
method used is applicable for mapping wind resources with a 
resolution of 0.5-10 km and the data needed are geostrophic wind, 
land and sea temperatures, topography, roughness and land-use. 
The simulations were made for 3 values of the geostrophic wind 


speed, for 16 wind direction sectors and for 4 representative 
months (January, April, July and October) to account for the 
seasonal variations in temperature. In addition, observations from 
six surface stations along with the National Center for Environ- 
mental Predictions (NCEP hereafter) reanalysis fields at the 
850 hPa pressure level were used to weight the influence of the 
geographical variation of the geostrophic wind into the results 
regarding the modeled wind climate [6]. 

Mederos et al. [7] presented a methodology to construct an 
offshore wind atlas for Canary Islands based on observations and 
model simulations. The authors produced mean 10-m wind speeds 
from daily wind data provided by QuickSCAT satellite for a 4 year 
period (2003-2006) at a resolution of 20 km and from 6-h winds 
obtained from HIRLAM model for the period 2002-2006 at 9-km 
resolution. The mean model values have been corrected based on 
the satellite data by establishing maximum and minimum values 
at 115% and 85% of the calculated mean value and values outside 
these limits were corrected and restricted to them. The “corrected” 
mean wind field of HIRLAM was used to “correct” with the same 
methodology the mean wind field produced from 3-h winds from 
MM5 model for the period 2004-2006 at 3-km resolution. Finally 
the resulting mean wind speed and mean direction at 10 m height 
with 3 km resolution are corrected with data from 40 anemometer 
stations. The authors concluded that a great offshore wind power 
resource is available in the Canary Islands. 

Gaston et al. [8] presented a wind resource map of Spain that 
was built using Eta/Skiron NWP model. The authors analyzed 365 
simulations covering one year (2006) at a horizontal resolution of 
0.1 x 0.1 deg. The presented wind atlas only focus on the statistical 
analysis of winds at 10-m. Following a similar methodology Costa 
et al. [9] developed a wind atlas for Portugal using model 
simulations with MM5 model for a whole year (1999) using four 
nested domains with 81, 27, 9 and 3 km grid resolution. For the 
model validation, observations from four surface meteorological 
stations were used. The results of the one year simulations were 
corrected with the mean deviation factor between the mean wind 
speed in 1999 and the long term mean wind speed recorded at 
4 anemometric stations, in order to consider the intra annual 
variability of the wind. The analysis focuses at three height levels 
(10-m, 60-m and 80-m). 

The wind atlas of Italy [10], first developed in 2002 and refined in 
2006 is based on the use of the mass-consistent WINDS (Wind field 
Interpolation by Non-Divergent Schemes) flow model developed at 
the University of Genoa. The authors have calculated the frequency 
distribution of wind speed and direction at 5000 m above ground 
level, provided by the European Centre for Medium-range Weather 
Forecast (ECMWF) reanalyses fields for the period 1990-1999. 48 
simulations (for 3 values of the wind speed and for 16 directions) 
have been produced, for a number of sub-areas, in order to cover the 
whole area of interest. The model was applied to obtain preliminary 
wind maps from geostrophic wind data, taking into account oro- 
graphy and terrain roughness. It should be noted that this methodol- 
ogy does not take into account local meteorological effects that 
are not related to the wind aloft (e.g. land and sea breezes). Every 
simulation has been weighted using the corresponding statistical 
frequency of the climatological wind aloft. Finally the simulated 
mean annual wind fields have been corrected through a “correction 
factor” obtained taking into account measured data from over 400 
stations scattered all over Italy. The model was applied to a grid of 
about 1 x 1 km. The refined wind atlas comprises the whole Italian 
land surfaces as well as the offshore areas up to 40 km from the 
coastline. 

Tammelin et al. [11] prepared the Finnish wind atlas by 
applying a series of models. The authors used initial and boundary 
conditions from the ERA-Interim data base and applied HIRLAM 
hydrostatic model at a resolution of 7.5 km over northern Europe 
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in order to provide initial and boundary conditions to the non- 
hydrostatic model AROME with 2.5 km horizontal resolution over 
Finland. The authors produced simulations for a period of 48 
months in total (4 Januaries, 4 Februaries, etc.) that were selected 
as representative of the period 1989-2007. Moreover they have 
also simulated the windiest and the calmest months of the same 
period. The selection of the representative months was based on 
the fulfillment of the following criteria: the Weibull distribution 
of the wind speed should be as representative as possible for the 
selected period and in the 12 wind direction sections while also 
the distribution of wind directions should be also as representative 
as possible. In addition to the NWP models, the diagnostic down- 
scaling method Wind Atlas Analysis and Application Programme 
(WAsP) [12] with 250 m resolution has been also applied only over 
the areas most favorable for wind power production. Further 
information on the Finnish wind atlas can be found at: http:// 
www.tuuliatlas.fi/en/index.html. 

Finally the Canadian Wind Atlas has been constructed using a three 
step procedure consisting of: (a) classifying the wind regimes such as 
to produce a climate state database based on NCEP reanalysis data, 
(b) producing simulations with MC2 mesoscale non-hydrostatic model 
at 5 km resolution, and (c) downscaling with a microscale modeling 
technique (MsMicro) which uses the wind energy statistics to deter- 
mine wind patterns at microscale [13]. The method is presented in 
detail in [14] and has been recently used to produce a high resolution 
wind atlas for Thailand [15]. A detailed review of the methods for 
wind resource estimation can also be found in [16]. 

For Greece, the wind atlas produced in 2006 [5] and updated in 
2008 has been based on the use of measured winds over the Greek 
territory and the use of a mass-consistent model (namely CRES- 
flowPot, see [17]). The methodology is based on the assumption 
that the wind flow at high altitudes is not influenced by the 
surface boundary layer, while close to the surface the combined 
effect of the topography and the boundary layer is sufficient to 
determine the wind speed and direction at any given point. The 
wind flow has been defined for 8 directions and normalized by the 
wind at the upper boundary. The normalized values of the study 
area have been then converted to actual wind speeds for each 
direction independently. Indeed, using time-averaging informa- 
tion (probability density function of the wind direction) that 
resulted from the measurements, wind speed has been calculated 
for each point as an average value. The characterization of the 
wind climatology over Greece has been based on the use of 200 
measuring stations deployed over Greece for various periods by 
the Centre for Renewable Energy Sources (CRES), the Hellenic 
National Meteorological Service (HNMS) and the National Obser- 
vatory of Athens (NOA). The available datasets included short-term 
(18 months) and long-term (10 years) measurements at various 
heights and periods. The datasets have been analyzed and corre- 
lated using the WindRose software tool developed at CRES (http:// 
www.windrose.gr) [5]. The results on the latest version of this 
wind atlas are provided on a mesh with a horizontal resolution of 
150 m and they are shown at the official web page of the Greek 
Regulatory Authority for Energy (http://www.rae.gr/geo/). 


3. The typical wind year 


The performance of wind turbines, much like solar systems, 
depends both on the long term-climatic-characteristics as well as 
the short term variations of weather conditions. Therefore, to 
predict the energy generated from such systems in the design 
stage, appropriate weather data is required. When the task at hand 
involves high resolution simulations the period that has to be 
simulated also needs to be carefully considered both in terms of 


representativity and in terms of efficiency/cost. In that respect 
running a multiyear simulation with a NWP model should provide 
adequately representative results otherwise it would be comple- 
tely inefficient. A possible way around this would to be to run the 
simulation for a single year with typical/representative weather 
conditions. Apart from the long term this would also take into 
account any possible seasonal variability as well as the short term 
fluctuations. 

Methodologies for generating weather data sets appropriate for 
this kind of studies, named Test Reference Years, Typical Meteor- 
ological Years (TMY) etc., have been in existence for many years, e.g. 
[18,19]. These data sets usually consist of 1 year of hourly data for 
the various meteorological parameters; temperature, humidity, solar 
radiation and wind are almost always included. These methodolo- 
gies can be divided into two major categories: (a) those generating 
synthetic data sets with statistical properties reproducing those of 
the long-term dataset, (b) those using sub-periods from the actual 
data set - usually months - with statistical properties best resem- 
bling those of the long-term data for that sub-period [20-23]. 

The existing TMY data sets of this kind were mostly intended 
for building energy and solar system simulations, thus the relative 
importance of wind in the selection procedure has been small 
and therefore they are inappropriate for a study of wind power 
potential [24]. 

For that reason, and in order to support the objective of the 
current study, a typical year based only on the wind flow over Greece 
and the surrounding maritime areas has been defined. For the 
definition of the typical year, the ECMWF meteorological reanalysis 
fields of 10-m wind speed and direction for the 20-year period 1989- 
2008 have been used. This ECMWF-Interim database provides data 
with ~80 km spatial resolution, at 6-h intervals. The selected area, 
for which the typical meteorological year was built, is bounded from 
18 E to 30 E and from 34 N to 42 N. 

The typical wind year consists of months selected from indivi- 
dual years, concatenated to form a complete year. As for the present 
application - the development of a wind-atlas - the parameters 
used from the ECMWF-Interim database are the 10-m wind speed 
and direction at 6h intervals for the aforementioned period. The 
methodology used in order to define the typical months is based on 
that introduced by Hall et al. [18] and proposed in [21] using the 
Filkenstein-Schafer statistic [25]. Indeed, the cumulative distribu- 
tion functions (CDFs) based on 160 wind bins (20 bins for the wind 
speed by 8 bins for the wind direction) have been calculated for 
each month, for each specific year of the 20 year period as well as 
for the entire 20 year period (long term CDF) for the same month. 
Then the wind CDF (speed and direction) for each month for each 
specific year is compared with the respective CDF of the long-term 
20 years composite. From the 20 months, the one for which the sum 
of the absolute differences between the two CDFs is minimal, is 
selected for the typical wind year. 

According to the aforementioned methodology the typical 
meteorological year constructed includes the monthly sequence 
given in Table 1. 


Table 1 
List of the 12 months and respective years that define the typical wind year. 


Month Year Month Year 

January 1999 July 2008 
February 1996 August 1997 
March 2000 September 2008 
April 1995 October 1994 
May 1998 November 2007 
June 2006 December 2003 
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4. Model simulations of the wind flow 
4.1. Modeling strategy 


As already presented in Section 3, model based wind atlases are 
based on the use of NWP models (such as HIRLAM, MC2, Eta, 
MM5, etc.) or microscale models or a combination of them. The 
microscale models are diagnostic models that are dedicated to the 
study of atmospheric processes that do not exceed a couple of 
kilometers and range from the simple versions (mass-consistent 
models, or more evolved linearized Navier-Stokes equation mod- 
els, e.g. WAsP, WINDS, MsMicro) to the advanced and computa- 
tionally demanding CFD models (e.g. MIUU) [16,26,27]. 

In the frame of this study the authors decided to use a NWP 
model for the construction of the wind atlas over Greece as the 
existing atlas has been based on the use of a mass-consistent 
model that is unable to include atmospheric parameters such the 
influence of atmospheric stability on the definition of the flow in 
the boundary layer, the development of local circulations (such as 
the sea breeze, katabatic winds, etc.) and the effect of Coriolis 
force. The advantage of the use of NWP model for the definition of 
the wind atlas is that the same model can be eventually used for 
the prediction of wind potential over an area (combined with a 
microscale model) creating thus a consistent strategy. 

The model based wind atlas of this study is based on the use of 
the NWP model MM5. MM5 model is a non-hydrostatic, primitive 
equation model using terrain-following coordinates [28]. Several 
physical parameterization schemes are available in the model 
for the boundary layer, the radiative transfer, the microphysics 
and the cumulus convection. For both the operational chain of 
MM5 model at NOA but also for the model simulations for the 
wind atlas the combination of Kain and Fritsch [29] parameteriza- 
tion scheme with the highly efficient and simplified microphysical 
scheme proposed by [30] have been used as they provide the most 
skillful forecasts of accumulated precipitation over Greece [31,32]. 
Concerning the choice of the Planetary Boundary Layer (PBL) 
scheme, the scheme proposed by Hong and Pan [33] is used as 
this has proved to have a consistently good behavior in the region 
{34]. This model has been verified for its forecasting skill over 
Greece by Kotroni and Lagouvardos [35]. 

Two nested grids have been defined and used for the model 
simulations (Fig. 1). Grid 1 has 10-km horizontal grid increment 
(230 x 210 grid points), covering the major part of Southestern 
Europe, and the Balkans. Grid 2 has 2-km horizontal grid incre- 
ment (431 x 436 grid points), covering the Greek territory and all 
the Greek islands and the adjacent water bodies. The horizontal 
extension of the defined grids is shown in Fig. 1a. In the vertical, 
forty unevenly spaced full sigma levels are selected. 

For each day of each month of the typical wind year, four 
simulations are performed starting at 0000, 0600, 1200 and 1800 
UTC with a 6-h duration for each simulation. All data from the 
6-h simulations sequences are stored for further analysis at 1-h 
intervals from t+1 up to t+6 h of simulation. The ECMWF gridded 
analysis fields and 6-h interval, at 0.5-degree lat/lon horizontal 
grid increment, are used to initialize the model and to nudge the 
boundaries of Grid 1 during the simulation period. 


4.2. Verification of the model simulations 


In order to evaluate the model skill to provide accurate wind 
forecasts, a verification procedure is undertaken for each month 
of the defined typical year. The observational data include six 
weather stations distributed across Greece (Fig. 1b). All stations 
are located at airports and operated by the Hellenic National 
Meteorological Service for periods longer than the examined one 
(1989-2008). The wind speed value at the model grid point closest 


to the observational site has been used for the verification, while 
the four model points surrounding the observing station were also 
considered. This procedure permits to reduce the possibility of 
penalizing the model performance due to inconsistencies between 
model and actual topography. 

The initial processing of the data, apart from consistency 
checking, included reduction of the model 10-m wind speed values, 
to the height of the observing stations, as well as binning of all data 
(modeled and observed) to 1 ms~! intervals. This treatment, even 
though deviating from the 0.5 m s~! recommended by the Inter- 
national Electrotechnical Commission (IEC) for wind turbines power 
performance measurements [36], was imposed by the fact that 
observed wind speeds were provided as integers. Furthermore, 
there are operational wind speed limits for power generation with 
wind turbines [37] that should be accounted for in the comparison. 
Thus, apart from the original, henceforth “complete” data set, we 
created a second dataset named “operational” by binning all data 
below or equal to a cut-in wind speed of 3 ms~! into a single bin. 
The cut-off speed (25 m s7!) was never reached in the analyzed 
dataset, so no specific treatment was necessary. 

In order to select the model point that best simulates the wind 
data observed in each meteorological station, contingency tables 
for each pair of observed and modeled data sets were created and, 
for each contingency table, the Gamma-statistic of association as 
well as the associated standard error were calculated [38]. Accord- 
ing to these statistics, wind speeds from any one of the 5 model 
points surrounding the observing station, represent the observed 
data equally well and in a statistically significant manner, either 
considering the complete or the operational data sets. 

Wind direction, even though not of interest in this particular 
context was also checked using contingency tables of observed 
and modeled data binned according to 9 (per 45°+ calms) or 5 (per 
90°+ calms) sectors. Only half of the 9-sector comparisons were 
statistically significant, whereas for 5-sectors only one station 
(Larissa) was not significant. However, as in the case of wind 
speed, there was no clearly discernible pattern as to which model 
points better represent the observations at each station. 

Based on the above, the model point closest to the observing 
station was used for the model performance evaluation. Basic 
model evaluation statistics for the model point closest to each 
station, and for both the complete and the operational data sets, 
are presented in Table 2. The model, on average over the entire 
wind speed range (complete dataset), overestimates wind speeds 
at mainland stations (Aktion, Alexandroupoli and Larissa), while it 
slightly underestimates wind speed at the island station of Limnos. 
When the operational data set is considered (wind speeds less or 
equal to 3 ms~! are removed) only slight overestimation remains 
at the mainland stations Alexandroupoli and Larissa, whereas 
wind speed is slightly underestimated at all other stations. The 
rest of the model performance statistics, also improve when we 
consider the operational data set instead of the original one, with a 
MAE around 2.2ms~—!, and Pearson correlation coefficients, r, 
exceeding 0.4 except for Larissa and Aktion. 

The typical point-to-point comparison presented above is 
useful in a general case, however it might not be the best choice 
in the specific context of our work. Since our goal is to provide 
credible information about the wind energy availability across 
Greece, it seems appropriate to focus on the comparison of the 
distributions describing the modeled and the observed wind 
speed data, rather than the data itself. The first logical step was 
to compare the empirical distributions of the modeled and the 
observed wind speeds. However empirical distributions are prone 
to reflecting the irregularities at the end points of the data, and 
in other places where the data is sparse, thus increasing the risk 
of a falsely negative ruling on the similarity of the modeled and 
the observed distributions. Furthermore, since our entire analysis 
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Fig. 1. (a) The two nested grids used for the simulations, (b) Inner model domain with the location of the surface stations used for the model verification: Alexandroupolis 
(A), Lemnos (L), Larissa (La), Aktion (Ak), Souda (S) and Heraklion (H). Model topography height (in m) is also shown. 
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Table 2 
Model evaluation statistics: Pearson correlation coefficient (r), Mean Bias Error 
(MBE), Mean Absolute Error (MAE) and Root Mean Square Error (RMSE) for the 
model point closest to each station, for the complete (C) and the operational 
(O) datasets. 

r MBE (m s7!) 


MAE (ms™!) RMSE (ms~!) 


(a oO E o Cc (0) C ie) 


Alexandroupoli 0.49 0.47 1.4 0.4 24 2.2 3.0 2.8 


Limnos 0.53 0.57 -0.7 -10 2.7 2.5 3.3 34) 
Larissa 0.18 015 3.0 05 34 2:1 4.0 2.8 
Aktion 0.33 0.21 15 -01 2.7 22 33 2.8 
Souda 0.41 0.43 0.2 -0.7 24 23 3.0 29 
Heraklion 0.31 041 -0.2 -10 24 24 3.0 31 


would be based on the fitted Weibull distributions for the model 
grid points, we preferred to use those as the basis of the compa- 
rison. An initial comparison using the Pearson correlation coeffi- 
cients, r, indicated very close agreement (r > 0.95 in all cases) 
between the CDF's of the distributions fitted to the observed and 
modeled wind speed. Moreover, the non-parametric Kolmogorov- 
Smirnov test for the equality of two probability distributions [39] 
was applied to the Weibull distributions fitted on the observed 
and the modeled data for each station and for both the complete 
and the operational data sets. 

Table 3 presents the resulting D statistics of the Kolmogorov- 
Smirnov test, as well as the respective critical values for the 5% 
level of significance. Values of D less than the critical level suggest 
that at the 5% level of significance, there is not enough evidence 
to reject the hypothesis that the two samples follow the same 
distribution. For the original set, at Alexandroupoli and Larissa the 
Kolmogorov-Smirnov test suggests that the modeled wind speeds 
do not follow the same distribution as the observed ones. However 
when the operational data set is examined (wind speeds between 
4 and 25m s`!) this discrepancy is removed and the Weibull 
distributions fitted on the modeled data are the same as those 
followed by the respective observed data. 


5. Statistical analysis of the wind 
5.1. Mean wind speed 


The Wind Atlas includes monthly wind statistics for seven height 
levels: 50 m, 75 m, 100 m, 125 m, 150 m, 200 m and 400 m, and for 
a total of ~ 190.000 points, each representing an area of 2 x 2 km’. 
The wind field has been analyzed at annual and monthly time 
scales, calculated from the hourly values. The area of the analysis 
covers Greece and the surrounding waters. In the following, the 
analysis focuses at the 50 m level winds. Fig. 2 presents the monthly 
average wind speed at 50m for one representative month 
per season. Indeed, during January the mean wind speed is less 
than 3 ms~! over the plains of continental Greece while over the 
medium and high altitudes the mean wind speed increases to more 
than 6 ms~! (Fig. 2a). Among the continental areas, Northeastern 
Greece, Attica and the adjacent island Evia, southern Peloponissos 
and the southern part of Crete present the higher mean wind 
speeds. Among the maritime areas, the Aegean Sea is much windier 
than the Ionian Sea, with mean wind speeds ranging from 7 ms~! 
up to more than 9ms7! locally. Increased wind speeds are also 
found in the entrance of the Corinthian Gulf, which is a narrow 
sea-level passage surrounded by a steep complex topography that 
consists of high mountains, elevated and sea-level gaps/straits. This 
is in agreement with the findings of Koletsis et al. [40] who found, 
in their wind climatology of the area that the most intense wind 
events occur during the winter season. 


Table 3 
Kolmogorov-Smirnov test, D statistic, the respective critical values for the 5% level 
of significance. 


D 

Ç (0) 

Derit=0.119 Derit=0.132 
Alexandroupoli 0.193 0.110 
Limnos 0.104 0.100 
Larissa 0.236 0.020 
Aktion 0.077 0.039 
Souda 0.067 0.067 
Heraklion 0.063 0.063 


During spring, (Fig. 2b) the mean wind speed over Greece 
presents much less variability than in winter, with winds of the 
order of 3-4 ms~! over the plains, increasing to 4-6 ms~! with 
height, while over the maritime areas the prevailing wind are of 
the order of 5 m s™!. 

During summer (Fig. 2c), while over land the mean wind speed 
pattern is quite similar to that of spring, over Attica, Crete and the 
maritime areas and the islands of the Aegean sea the wind speed 
is considerably increased. Indeed, during summer northern sector 
winds are blowing over the Aegean Sea, the Etesian winds, which are 
mainly north-easterly in the northern Aegean, northerly in the 
central and southern Aegean and tend to become north-westerly 
near the southwestern Turkish coasts. The air masses regularly 
originate from the region of southern Russia and the Caspian Sea 
and they are dry and relatively cool, contributing to the decrease of 
surface temperature and the moderation of summer heat and 
discomfort [41,42]. The sustained wind speed associated with the 
Etesians often attain near-surface values exceeding 15 ms~! [42]. As 
shown in Fig. 2c, the mean wind speed at 50-m, over the Aegean, 
increases to 7-8 ms~! and locally reaches more than 11 ms~!. 

Autumn is the rainiest season in Greece which is affected by the 
passage of low pressure systems associated with southerly winds 
[43]. The mean wind speed distribution at 50-m is relatively close to 
that during winter with the windiest areas found over the Aegean, 
Attica, Evia, the Corinthian Gulf area and Thrace in northeastern 
Greece (Fig. 2d). 


5.2. Wind potential 


The model simulated winds at hourly intervals at 50-m height 
have been used in order to determine the wind potential. If fv) is 
the probability density function (PDF) of wind speed v, the wind 
potential P can be expressed as:. 


1 oc 
Pu) = 5p f vf (v)dv (1) 


where p is the air density. This integral can be evaluated over a 
range of wind speeds and in the frame of this study it has been 
evaluated for the range of operation of wind turbines that usually 
ranges from a cut-in speed of 4m s7! up to the cut-off speed of 
25m s™' witha step of 0.2 m s™!. 

The PDF of wind speed is important in numerous wind energy 
applications. [44] presented an extensive review of the available 
methods for its estimation. Among a variety of PDFs, the two 
parameter Weibull probability density function is the most widely 
used and accepted in the specialised literature on wind energy 
and other renewable energy sources [45]. Indeed, the wind speed 
probability density function is given by 
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Fig. 2. Monthly mean wind speed (in m s7!) at 50-m for (a) January, (b) April, (c) July and (d) October. 


where c (in m s7!) is the Weibull scale parameter and k (dimen- 
sionless) is the Weibull shape parameter that characterizes the 
asymmetry of the probability function. There are various methods 
to determine these two parameters. 

In the frame of the present work, for the estimation of the two 
parameters the least square method has been used (see [44]). 
Indeed, the shape parameter k can be estimated through the 
following equation: 


N N N N N 2 
k=fn $ vx- 2x) (r) n Zx- (x) 
t= i=1 i=1 i= i= 


where 
Y; = In{ —In(1 —P;)} and X; = In(Vmax.i) 


(3) 


and P; are the components of the vector P which contains the 
cumulative relative frequencies of the N intervals in which the n 


wind speeds of the sample have been classified, Vmaxi are the 
maximum wind speed values within each of the N intervals. The 
scale parameter c is estimated from the following equation: 


b 
c=exp |- a (4) 


where b is estimated following: 


N N 2 N N N 2 N 2 7 
ba { 3 v3 x8 (3x1) ($ xen) hn 3 -h E x) 


(5) 


Based on the above, once the wind potential is estimated a 
classification was made following [46]. Indeed, six classes of wind 
potential have been used as presented in Table 4. 

Fig. 3 presents the mean yearly wind potential at 50-m that has 
been evaluated from a cut-in speed of 4 ms~! to the cut-off speed 
of 25 ms! with a step of 0.2 ms~! based on Eq. (1). It is clear 
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Table 4 
Classification of wind potential density (adopted from the Texas State Energy 
Conservation Office). 

Wind potential density (W m~?) 


Wind Power Class Commercial viability 


1 0-200 Very poor 
2 200-300 Poor 

3 300-400 Marginal 
4 400-500 Good 

5 500-600 Very good 
6 > 600 Excellent 
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Fig. 3. Wind potential (in W m~?) at 50-m. 


that in the largest part of the Aegean sea the wind potential is 
larger than 500 W m~”, reaching locally 700 Wm~?. Over the 
Ionian the wind potential is of the order of ~300 W m~°. Locally 
the wind potential exceeds 400 W m~? over the mountain ranges 
as well as over the Corinthian Gulf, Thrace, Attica and Evia, 
southeastern Peloponnisos and Crete. Accordingly these areas are 
characterized by wind power class 3 or higher as depicted in Fig. 4. 
A zoom-in over the areas with the highest wind classes, namely 
Crete (Fig. 5a), and Attica and Evia (Fig. 5b) is given. It is evident 
that eastern Crete has a higher class wind power that locally is 
excellent (class 6), while southern Crete presents a better wind 
power class compared to the northern areas. This is related to the 
complex topography of Crete island that produces gap winds and 
acceleration of the wind flow over southern Crete [47,48]. The 
wind power class exceeds 5 over a large part of Southern Evia 
island, while offshore southern Evia the wind power class rates as 
excellent. 

With the aim to investigate the effect of the topography height 
and slope on the wind distribution and consequently on the wind 
potential density, the cumulative frequency of wind power class 
for 7 classes of topography height and 5 classes of slope have been 
calculated. Namely, the selected topography bins are: < 200, 200- 
400, 400-600, 600-800, 800-1000, 1000-1500, > 1500 m) while 
the slope bins are: 0, 0.1-2, 2-5, 5-10, > 10). At this point it 
should be noted that the topography heights used are those used 
by the model at the 2 km x 2 km grid increment and thus they are 
smoothed compared to the actual topography heights. Moreover 
the sea points have been analyzed separately. The analysis 
(Table 5) shows that: 
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Fig. 4. Wind power class at 50-m. 


- Almost 80-85% of all the land points with topography heights 
up to 1500 m are characterized by a wind power class of 1 
(< 200 W m~?) while the percentages of class 2 range from 8- 
13%, for class 3 range from 2.9-4.5% and only about 0.8% of the 
land grid points are characterized by a wind power class of 6. 

- Grid points with heights above 1500 m present a higher wind 
potential density. Indeed ~60% of the points of this height 
bin are characterized by a wind power class of 1, while the 
percentages for classes 2 and 3 are 30% and 6%, respectively. 

- The dependence of the wind potential density on the slope is 
more evident. The grid points of O slope are characterized by 
very poor wind potential density at a percentage of ~90%, 6% 
are classified to the class 2, ~2% to class 3 and only 1% of the 
points are classified in class 4 or higher. The grid points with 
slope up to 5%, present a better commercial viability. About 
11-15% of them are in class 2, and almost 6% in class 3. Finally 
for the grid points with slopes larger than 5% the wind poten- 
tial density distribution is even better. Indeed ~15% of the 
grid points are characterized by class 2 wind power, ~6% by 
class 3 and ~6% of the points are characterized by a wind class 
category of 4 or better. 

- Analysis of the sea points showed that the offshore wind 
potential over the Greek seas is quite important as almost 
70% of the marine grid points are characterized by a wind class 
of 3 or better, while ~ 34% of the grid points are classified in the 
wind power class of 4 or better. 


6. Concluding remarks 


In the frame of the present paper a Wind Atlas over Greece has 
been developed based on simulations performed with the numer- 
ical weather prediction model MM5 that has been actively utilized 
at the National Observatory of Athens since 2000. The atlas 
contains 12 months of weather model simulations where the 
winds and other variables are calculated for grid boxes 2 x 2 km? 
in size. 

The methodology used in order to select typical months to 
form a typical wind year was based on the Filkenstein-Schafer 
statistical method [25] and it is described with details in Section 2. 
The typical wind year was simulated by MM5 in 6-h sequences, 
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Fig. 5. Wind power class at 50-m. Zoom in Fig. 4 for (a) Crete island, (b) Attica and Evia areas. Only areas with class 3 or higher are shown. 


storing the atmospheric state at 1-h intervals from t+1 up to 
t+6h of simulation. The initial state of each 6-h sub-period and 
the required boundary conditions were provided by the European 
Centre for Medium-Range Weather Forecasts (ECMWF) analysis. 
The results of the model simulations have been verified against a 
number of surface meteorological stations while also annually and 
monthly statistics at various heights are produced over Greece and 
the surrounding maritime area. 


The analysis showed that among the continental areas, North- 
eastern Greece, Attica and the adjacent island Evia, southern 
Peloponissos and the southern part of Crete present the higher 
mean wind speeds, while over the maritime areas, the Aegean Sea 
is much windier than the Ionian Sea. As for the seasonal variability 
of the mean wind speed, it was evident that during summer the 
wind speed is considerably increased over Attica, Crete and the 
surrounding maritime areas and the islands of the Aegean sea. 
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Table 5 
Distribution of wind power classes (in %), in the studied domain according to 
5 slope and 7 terrain height categories, as well as for the sea points. 


Wind power class 


al 2 3 4 5 6 


Slope 0 90.9 5.9 22 1.0 0.0 0.0 
0.1-2 80.7 11.9 4.3 2.0 0.7 0.4 
2-5 78.6 11.8 5.5 2.3 Wl 0.7 
5-10 72.8 15.1 5.8 3.2 1.4 17 
>10 73.6 13.9 6.1 2.9 13 2.2 
Terrain Height < 200 85.7 10.0 3.2 0.8 0.2 0.1 
200-400 81.3 112 45 1.8 0.7 0.6 
400-600 80.4 11.1 4.5 2.4 0.9 0.8 
600-800 83.8 8.2 3.7 23 11 0.8 


800-1000 85.2 8.4 29 1.8 1.0 0.8 
1000-1500 80.9 12.9 4.0 13 0.4 0.6 
> 1500 62.1 29.2 5.9) 13 0.6 0.9 
SEA 8.3 223 35.8 21.6 9.3 2.7 


Moreover, in the largest part of the Aegean sea the wind 
potential is larger than 500 W m~?, reaching locally 700 W m~?. 
Over the Ionian the wind potential is of the order of ~300 W m~2. 
Locally the wind potential exceeds 400 W m~? over the mountain 
ranges as well as over the Corinthian Gulf, Thrace, Attica and 
Evia, southeastern Peloponnisos and Crete. These latter areas are 
characterized by wind power class 3 or higher. 

Finally the investigation of the effect of the topography height 
and slope on the wind distribution and consequently on the wind 
potential density showed that the latter is more dependent on the 
slope than on the topography height. The analysis of the sea points 
showed that the offshore wind potential over the Greek maritime 
areas is quite important as almost 70% of the marine grid points 
are characterized by a wind class of 3 or better. 

At this point it should be mentioned that for wind resource 
estimation and siting a detailed description of the terrain, namely 
the surface roughness, orography and obstacles, is required. In that 
case the Weibull distributions calculated over Greece could be 
used as input to a microscale model such as WaSP [12,49] together 
with detailed description of the terrain topography in order to 
estimate local wind climate for wind energy applications. Thus, 
it is in the authors plans to build a web-based interface that will 
make available the current wind atlas results for use for siting 
purposes. 
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